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Abstract 

The relationship between the Coulomb displacement energy for the A=8, 
J=2+, T=l state and the low-energy astrophysical factor for the 
^Be(p,7)^B reaction is discussed. The displacement energy is interpreted in a 
particle-hole model. The dependence of the particle displacement energy on 
the potential well geometry is investigated and is used to relate the particle 
displacement energy to the rms radius and the asymptotic normalization of 
the valence proton wave function in '^B. The asymptotic normalization is used 
to calculate the astrophysical S\-j factor for the ''Be(p,7) reaction. The rela- 
tionship to the ^Li(n,7) reaction, the ^B quadrupole moment, radial density, 
and break-up momentum distribution are also discussed. 



1 



Typeset using REVTJj^X 



I. INTRODUCTION 



Measurements of high energy solar neutrinos in the Homestake ^ and Kamiokande II 
and III ^ experiments have found a significantly smaller number of solar neutrinos compared 
to that expected from the standard solar model. ^ The /?+ decay is the main source of 
these these high energy solar neutrinos. is formed by the ^Be(p,7)^B reaction at a center 
of mass energy of about 20 keV. Its cross section is conventionally expressed in terms of the 
5*17 factor. ^ Sn is known only from an extrapolation of data at energies above 100 keV. Also 
the highest precision data disagree, with those of Parker ^ and Kavanagh et al. ^ being about 
30 percent higher than those of Filippone et al. ® and Vaughn ct al.. ^ This has lead to the 
investigation of other ways to determine the Si'j factor such as the Coulomb dissociation, ^ as 
well as other properties of ^B which are indirectly related to the Sn factor such as its total 
reaction cross section, ^ quadrupole moment, and break-up momentum distribution. 

Much of the theoretical work has been based upon using a potential model to gener- 
ate the single-particle wave function for the most loosely bound proton, and then combining 
this with a shell-model calculation of the spectroscopic factor to obtain the asymptotic nor- 
malization of the wave function. This current work was initiated by the observation that 
the potential model parameters (e.g. the radius R and diffuseness a of the Woods-Saxon 
potential) used in previous calculations are taken from some "standard" sets based upon 
nucleon-nucleus scattering optical potential analyses, and that these "standard" parame- 
ters were not obviously appropriate to the particular case of ^Be plus protons. We thus 
investigated the extent to which the potential parameters could be determined from the dis- 
placement energy of the A=8, T=l state as given by the binding energy difference 
between the ^B and ^Li ground states. In addition, we also examined several quantities of 
related interest including the ^Li(n,7) cross section, the Q moment, the density distributions, 
and the momentum distributions. 
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II. COULOMB DISPLACEMENT ENERGY 



The relationship between the root-mean-square (rms) radius of the valence proton 
and the displacement energy is well known. ^^'^^ Qualitatively, the larger the rms radius the 
smaller the displacement energy since the valence proton is further from the core protons. 
This leads to the Thomas-Ehrman effect in which loosely bound valence nucleons have a 
relatively smaller displacement energy compared to more tightly bound nucleons. In addi- 
tion, for a fixed binding energy, nucleons in a low £ state have a smaller displacement energy 
compared to those in a high I state because the centrifugal barrier for high I results in a 
smaller rms radius. For heavier nuclei (A>16) with a relatively simple shell-model config- 
uration, many quantitative calculations of the displacement energy have been carried out. 
At the beginning of the Odls shell (A=17, T=l/2) and the Oflp shell (A=41, T=l/2), the 
observed displacement energy is found to be about 10 percent larger than that obtained 
from calculations which involve only the lowest shell-model configuration together with the 
Coulomb interaction. These cases include those for high I values whose Coulomb shift is not 
very sensitive to the potential-well geometry. This "Nolen-Schiffer" (NS) anomaly appears 
not to have a simple explanation but is due to a combination of core-polarization, high-order 
configuration mixing, and charge-asymmetric strong interactions effects. The systematics of 
the displacement energies of light nuclei can also be semi-quantitatively accounted for with 
a simple potential-well geometry. 

In this paper we apply what has been learned about the displacement energy system- 
atics to constrain the shape of the potential for the ^Be plus valence proton system. The 
A=8, T=l displacement energies are more complicated than the T=l/2 systems usually 
considered, since at least two nucleons are involved. However, the Op shell-model wave func- 
tions predict a relatively simple structure for A=8, T=l. The spectroscopic factors (C^S) 
for (A=8, 2+, T=l) ^ (A=7, 3/2", T=l/2) and (A=9, 3/2", T=l/2) ^ (A=8, 2+, T=l) 
are experimentally and theoretically both near unity. ^^'^^ Thus, can be considered as 
a proton-particle neutron-hole configuration relative to ^Be. The mirror nucleus ^Li is a 
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proton-hole neutron-particle configuration. ^Be is itself "deformed" and thus exits in both 
the ground state and 2+ first excited state relative to the particle-hole configuration. 

The displacement energy of the particle-hole state is simply the sum of the particle 
and hole displacement energies (the Coulomb particle- hole interaction is zero): 

As = Aft + Ap, (1) 

where Ag is the A=8, T=l displacement energy (the binding energy difference between ^Li 
and *B). In the simplest approximation A^ is the A=7, 3/2~, T=l/2 displacement energy and 
Ap is the A=9, 3/2", T—1/2 displacement energy. The experimental values are A7=1.645 
MeV and A9=1.851 MeV which gives Ay -|- A9=3.496 MeV compared to the experimental 
A=8 value of A8=3.540 MeV. So the simplest model works rather well. There are however 
several reasons why this is not exact. One is that the actual shell- model configuration is 
a little more complicated than just particle-hole. Also the Thomas- Ehrman effects may be 
significant since the proton is unbound by by 186 keV in and bound by 138 keV in ^B. 

Since we are particularly interested in the properties of the valence proton in ^B we 
will proceed as follows. First Eq. (1) will be modified to read 

As = Aft + a; + A,„, (2) 

where A^^ will take into account the Op shell-model structure beyond particle-hole. Since 
the hole state is relatively tightly bound, the Thomas-Ehrman shift will be small and we 
will take Aft = Ay = 1.645 MeV. A^ indicates the displacement energy for the particle in ^B 
which may differ from Ap = Ag because of the Thomas-Ehrman shift. 

The shell-model correction Agm was calculated by using the Coulomb plus charge 
asymmetric interaction of Ormand and Brown within a full Op shell-model basis. The 
matrix elements were calculated with harmonic-oscillator radial wave functions. The results 
for the displacement energies are: Aft = A7 = 1.719 MeV, As = 3.656 MeV and Ap = 
Ap — Ag — 1.873 MeV. This gives A^^ = 64 keV. Since harmonic-oscillator radial wave 
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functions are used the shell-model calculation should not necessarily be in good agreement 
with experiment, however, it should provide an estimate for the A^^ correction. Thus the 
displacement energy of the A=8 particle state is 



a; = As - A;, - A,^ = 



1.831 MeV. 



(3) 



We are interested on how A^ 



depends upon the potential geometry, and on how the 



potential geometry effects the rms radius and the astrophysical 5*17 factor. To investigate 
this we will calculate the direct part of the Coulomb shift in a Woods-Saxon geometry. This 
means that we calculate the single-particle binding energy with and without the one-body 
Coulomb potential and take the difference to obtain the Coulomb displacement energy. The 
Woods-Saxons potential has the usual form of central plus spin-orbit plus Coulomb terms. 
The central potential has the form 



where R is the radius, a is the diffuseness. We use the reduced mass in the kinetic energy 
operator. The spin-orbit potential is the usual derivative form with the same geometry 
as the central (we will show below that the spin-orbit potential is not important for our 
analysis). The Coulomb potential is obtained from the density distribution of four tightly- 
bound protons (two in Os and two in Op) obtained with another Woods-Saxon potential 
which is constrained to reproduce the experimental rms charge radius of 2.52 fm for ^Be. 

We thus investigate the dependence of the Coulomb energy, the rms radius and the 
astrophysical Sn factor on the radius R and the diffuseness a of the central potential. In 
all cases we vary R and a and fix Vg in order to reproduce the proton separation energy of 
138 keV. We calculate the direct Coulomb displacement energy Ao and relate it to A^ by 
making several additional corrections, all of which are relatively small: 



V{r) = Vo{l + exp[r-R\/a\} 



-1 



(4) 




+ ^ex + ^so + ^vp + ^np + ^NS- 



(5) 
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Agj. is the exchange correction. A value of Aea;=-125 keV was obtained from the harmonic- 
oscillator shell-model calculation for Ag by comparing the results with and without the 
Coulomb exchange terms. The next three terms are the relativistic spin-orbit correction 
(Eq. 21 in Ref 12), the vacuum polarization correction (Eq. 4.23 in Ref 13) and the proton- 
neutron mass difference correction (Eq. 4.29 in Ref 13), respectively. Our estimates for these 
are —30 kev, 12 keV and 14 keV, respectively. The Nolen-Schiffer correction, Ajv5 is perhaps 
the most uncertain. For A—YI and A=41, A^vs/Aq is about 0.10 and we see no reason why 
it should differ much from this for the beginning of the Op shell. Since Aq is about 1.8 MeV 
we will take Ajv5=180 keV. Thus we arrive at the empirical value of 



with which we will compare our calculations for the direct Coulomb shift of the valence 
particle. 

We show in Table I the results for Ao, the rms proton radius and the density of the 
valence proton at r = 10 fm, p(lOfm), as a function of R and a. Beyond the influence of the 
strong interaction (about 6 fm) the shape of the radial wave function is entirely and uniquely 
determined by the Coulomb plus centrifugal potentials. The only quantity which depends 
on the potential is the asymptotic normalization as represented, for example, by the value 
of p at r = 10 fm. Our p(r) is defined by the normalization: 



In order to show the correlation between Aq, rms and p(lOfm) we plot values obtained for 
these quantities in pairs in Figs. 1 (Aq vs rms), 2 [Ag vs p(lOfm)] and 3 [rms vs p(lOfm)]. 
We find that there is strong correlation between them, which implies that from a knowledge 
of any one of them (in particular Aq) we can infer a range of values for the other two. 

The results in Table I were obtained for a Op3/2 valence particle. For R—2A fm we 
show the results with and without the spin-orbit potential, and from this comparison it is 



Ao(empirical) = 1.780 MeV, 



(6) 




(7) 
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clear that spin-orbit potential is not important (as long as is fixed from the 138 keV 
separation energy). Thus to a good approximation, our results apply to both Op3/2 and 
Opi/2- 

The results shown in Figs. 1—3 show that there is a rather narrow band of points 
which relate the quantities of interest in the framework of the Woods-Saxon potential model. 
From Fig. 1 one finds that the value of Ao=1.780 MeV corresponds to a narrow range of rms 
valence radii from 4.5 to 4.7 fm, and from Fig. 2 one finds that the same Aq corresponds to 
p(lOfm) values in the range 7.2 to 8.0 fm~^. The correlation between the rms radius and p 
shown in Fig. 3 is even more parameter independent. 

The potential parameter values required fall into three sets with (i?, a) values of 
about (2.8, 0.51 fm) (set A), (2.4, 0.69 fm) (set B), and (2.0, 0.81 fm) (set C). Or analysis 
is consistent with any of these sets or any interpolation between them. These are compared 
in Table II to the parameters used in a number of other potential model calculations for 
5'i7. The results for Ao, rms and p(lOfm) for the potentials in Table II are shown by the 
labeled crosses in Figs. 1-3. The results with Tombrello parameters (T) are close to ours. 
They were obtained from the optical model analysis of 180 MeV proton scattering on Li and 
Be by Johansson et al.. The Johansson parameters were subsequently used by Aurdal^'^ 
and Robertson. Parameter set I from Barker is also close to ours. But his modifications 
of the parameters needed to reproduce the ^Li(n,7) cross section (see the discussion below) 
are clearly outside of our range. 

III. RESULTS FOR THE ASTROPHYSICAL Sn FACTOR 

The astrophysical Syj factor is related to p(lOfm) in the following way. At very low 
energies (< 50 keV) the ^Be(p,7)^B capture cross section is dominated by the El transition 
between the ^Be-|-p scattering states and the *B ground state, and it is almost exclusively 
determined by contributions coming form the external part of the scattering and bound state 
wave functions. The asymptotic behavior of the scattering states is uniquely defined (the 
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phase shifts being the hard sphere phase shifts - practically zero), and the asymptotic part 
of the wave function, which describes the ^Be+p relative motion in the bound state ^B, is 
proportional to the fixed Whittaker function. 



i^ssir) = ci 



r 



oo. 



(8) 



r 



where / = 1, 2 is the channel spin, r is the radial distance between '^Be and p, and k is the 
wave number corresponding to the ^B binding energy relative to the ^Be+p threshold, cj 
are the constants which are required to normalize the Whittaker function to the asymptotic 
^B wave function. Thus, at low energies the astrophysical 5'i7 factor of the capture reaction. 



depends only on c. ^^'^^ (here rj — e^ZiZ^/hu with Zi — 1, and Z^ — A is the Sommerfeld 
parameter). From the hard sphere scattering states one obtains: 



(The uncertainty coming mainly from the fact that the nucleon mass is not well defined in 
nonrelativistic quantum mechanics, is 1—2 percent. The value of the S factor at zero energy 
is roughly 0.4 eV-barn higher than at 20 keV.) Using this formula we can express Sn (in 
units of eV-barns) in terms of the valence proton density at any given asymptotic radius, 
e.g. at 10 fm. 



5i7(20kev) = 2.99-10^3/2(10 fm) S3/2 [(0:1,3/2 + 7ai,i/2)' + (a2,3/2 + 7a2,i/2)'], (H) 



were the a coefficients are determined from the transformation between the ^Be( J^) + £j — 
^B( J/) coupling and the channel-spin coupling: 



Si,{E) = a{E)Eexp[27iri{E)], 



(9) 



5i7(20keV) =36.5(c? + c^). 



(10) 



[Ji ij]'^ = E ^lAiJi ^ 1/2)' «) iV^, 



(12) 



I 
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and 7 is given by the ratio: 



_ [^1/2 ^i/2(10fm)] 

^ [^3/2 V'3/2(l0fm)]" ^ ^ 

In Eq. (13) 9j is the n = 0, i=l spectroscopic ampUtude and in Eqs. (11) Sj is the spectro- 
scopic factor. The amphtudes 9j are given by the reduced matrix elements of the creation 
operator, a'^: 

OM^ Jf) :j^j=^ , (14) 

When 9j is given without its Jj, Jf arguments as in Eq. (13), it corresponds to the ground 
state (Ji = 2) to ''Be ground state (J/ — 3/2) value. The spectroscopic factors take into 
account the additional center of mass correction factor, [Ai/{Ai — 1)] = 8/7: 



S> = (15) 
In Eq. (13), t/'(10fm) is the radial amplitude at r = 10 fm: 



Pj(10fm) =-01(10 fm). (16) 

In our case where £=1, Ji — 3/2, Jf — 2 and / = 1,2, the transformation coefficients 
are ai^3/2—a2,3/2—'^2,i/2—i/V^ and ai^i/2—-l/V^- In addition, to a good approximation, 
t/'i/2(10fm) = t/'3/2(10fm) = -0(10 fm). Hence Eq. (11) simplifies to: 

,5i7(20kev) = 2.99-10^(10 fm) S, (17) 

where S ^ Ss/2 + S1/2. 

In Table III we give the values of 9j obtained from Op shell model calculations with 
a variety of interactions which are appropriate for the lower part of the Op shell. These 
amplitudes are quite stable with respect to a reasonable range of interactions. 
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In Table IV we compare the theoretical spectroscopic factors with those obtained 
from reaction data for states of ^Li as well as those extracted from the observed widths 
of unbound states in ^Li and ^B. The spectroscopic factors for the unbound states are 
obtained from F^xp — SFsp where Fsp is the single-particle width for a resonance at the 
experimental separation energy in the potential geometry B (i? = 2.4 fm and a = 0.69 
fm). The decay data are observed to be in excellent agreement with theory. The reaction 
spectroscopic factor for the 2+ state is low with respect to theory, however, it depends upon 
the potential parameters used for the DWBA calculations (we have not attempted to repeat 
the DWBA calculations). The reaction spectroscopic factor for the 1"*" state is also low with 
respect to the decay spectroscopic factor of the mirror state. The data in Table IV support 
the Op shell-model calculation for S and the present potential parameters. 

Combining with our results of 5* = 1.15 ± 0.05 and p(lOfm) = (7.7 ± 0.4)10"^ fm~^ 
we obtain ,Si7(20keV) = 26.5 ± 2.0 eV-barns. 

IV. RELATIONSHIP TO THE 7LI(N,7)8LI CROSS SECTION 

Barker^® has pointed out that standard potential models tend to overestimate the 
experimentally well-determined low-energy ^Li(n,7)^Li cross section. He argued, that one 
should modify either the potential parameters or the spectroscopic factor to get agreement 
with experiment. To study this issue, we performed calculations for the ^Li(n,7)®Li reaction. 

Because there is no Coulomb barrier in this reaction, the inner parts of the wave 
functions have the same importance as the asymptotic parts. Thus the cross section does 
not depend solely on the asymptotic normalization of the bound state wave function. For 
the ^Li bound state we used the same potential parameters as for ^B, except a change in 
the potential depth to get the exact neutron separation energy of 2.033 MeV. For the "^Li+n 
scattering states, we modify the potentials to reproduce the experimental scattering lengths 
of the I — 1, and 1 — 2 channel spin states, respectively. 

The thermal ^Li(n,7)^Li cross section in our model, using potentials A, B, and C, are 
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78, 80 and 83 mb, respectively. The experimental thermal cross section is 45.4±3.0 mb. '^'^ 
Thus we also obtain an overestimation of the ^Li(n,7)®Li cross section as did Barker. He 
concluded that either the potential parameters or the spectroscopic factor has to be changed 
in order to agree with the experiment, and that these changes would bring the ^Be(p,7)^B 
5'i7(20keV) factor down to 16-17 eV-barns. The modifications to the potential are very large 
(R is changed from 2.39 to 0.53 fm or a is change from 0.65 fm to 0.27 fm) and we can see from 
Figs. 1—2 that these large changes are inconsistent with the Coulomb displacement energy. 
Thus we can exclude the possibility of radically changing the potential parameters. In the 
spirit of Ref 28, the only remaining possibility would be the reduction of the spectroscopic 
factor to about 0.71. But given the general agreement we obtain for the decay widths in 
Table IV, such a large change in the spectroscopic factor seems unreasonable. In fact such 
a drastic change would question the adequacy of the potential model itself. 

We would hke to point out that the discrepancy in the '''Li(n,7)^Li cross section could 
be resolved in a way which does not affect the '''Be(p,7)®B cross section. As mentioned, 
contrary to ^Be(p,7)^B, the inner part of the wave functions are important in the case of 
^Li(n,7)^Li. Although the reproduction of the scattering lengths fixes the external part of the 
scattering wave functions, the internal, off-shell, part is not well-constrained. For instance, 
if the inner node of the wave function were somewhat further outside than in the potential 
model, this would bring the '^Li(n,7)^Li cross section down. To illustrate that the node 
position is not well-defined, we show in Fig. 4 the inner part of the 1 — 2 scattering wave 
function of the standard potential (BI) of Barker (solid line) together with the scattering 
state obtained from the cluster model of Ref 29 (dashed hne) at Ecm—W keV. This change 
in the off-shell behavior is enough to reduce the '''Li(n,7)^Li cross section considerably. In 
fact, the dashed line of Fig. 4, together with the bound state of the standard potential of 
Barker results (after a, 1/v extrapolation from 10 keV) in a thermal cross section of 46.3 
mb, which is close to the experimental value. We emphasize again, that this modification in 
the off-shell behavior of the scattering wave functions has no effect on the ^Be(p,7)*B cross 



11 



section. 



V. RELATION TO THE QUADRUPOLE MOMENT 

The quadrupole moment, Q, of is related to the above calculations in the following 
way. Q is proportional to the matrix element of the E2 — r^y^^^ one-body operator whose 
reduce matrix element is given by a summation over products of many-body matrix elements 
times single-particle matrix elements : 

(J. II E2 II J,) = E iAM^k^^^flMt^,,,^ II E2 II f,Q, (18) 
o ,/ f V -I- i 

where \—2. The tz indicates a sum over protons and neutrons. By inserting a complete set 
of states (J/) of the A = 1 system between the a"*" and a, one can rewrite this as a sum 
over all Op shell states of the A=7 system: 

(J, II E2 II J,) = (2 J, + 1) E i-^y'^''OjM, Jf)Oj',tM: Jf) 

j,j',tz,Jf 

X I ^ i 0',iz II E2 II (19) 

J f Jf 

where the states Jf are in ^Be for i^^^proton and are in ^B for t^^neutron. The sum over 
Jf can be broken down into the term coming from the ground state of ^Be (referred to as 
the valence term, vp) and all other terms coming from excited states in ^Be and all states 
in ^B (referred to as Op core proton and neutron terms, pep and pen, respectively). The 
single-particle matrix elements are given by a geometrical term times the single-particle 
mean-square radius, (r^).^° The Q moment can thus be expressed in the form 



QCB) = -0.80 ^3/2^1/2 (r')^pep + 0.203 {r')pcpep + 0.183 {r')pcnen 

= 0.187(r2)„pep + 0.203 {r^)pcpep + 0.183 {r^)pcnen, (20) 
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where the CKI interaction was used for the spectroscopic factors. The numerical coefficient 
—0.80 is purely geometrical. The effective charges Cp and e„ take into account the non- 
Op parts of the wave functions which include Os to Op and Op to Odls proton excitations. 
For the remaining discussion we will use values of Cp — 1.35e and e„ = 0.35e, ^° although 
we realize that these are approximate values and that they may depend upon the binding 
energy. Our calculations for the radial matrix elements (with potential B) give (r^)^p = 
21.6 fm^, (r^)pcp = 8-1 fm^, and (r^)pcn = 7.8 fm^, and hence we obtain Qi^B) — 8.2 efm^, 
in reasonable agreement with the experimental value of 6.83±0.21 efm^. The Q moment 
for ^Li can be obtained by interchanging the labels for p and n in Eq. (20), and using our 
value of {r'^)vn — 14.0 fm^, to obtain Qi^lA) — 3.47 efm^ which is close to the experimental 
value of 3. 2 7± 0.06 efm^. The effective charges we use are approximate and it would very 
difficult to estimate the them more quantitatively. Even though our calculationed values 
are close to the experimental values, it is clear that the relationship to the valence rms 
radius and asymptotic normalization is quite complicated. In addition, we note that the Q 
moment depends upon the interference term ^3/2^1/2 whereas the tail density is determined 
by the combination ^3^2 + ^1/2 which is dominated by ^3/2. The term in Eq. (19) which is 
proportional to ^3/2^3/2 is zero because the 6j symbol vanishes (physically it is related to the 
vanishing of the Q moments in the middle of a single j shell) . 

VI. RADIAL DENSITIES AND MOMENTUM DISTRIBUTIONS 

The radial densities obtained in our potential B calculations are shown in various 
ways in Figs. 5-7. For the sake of simplifiction it will be assumed in this section that there is 
a single valence proton with a binding energy of 0.138 MeV, four tightly bound core protons 
(two in the Op shell with a separation energy of about 6 MeV and two in the Os shell with a 
separation energy of about 16 MeV) and three tightly bound neutrons (one on the Op shell 
with a separation energy of about 8 MeV and two in the Os shell with a separation energy of 
about 18 MeV). The actual situation is a httle more complicated than this because S — 1.15 
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and because there is also some parentage of the protons in to the first excited state in 
^Be which will result in some leakage of the core protons to larger radii. The results here 
are thus more qualitative than those given above for Sn. 

The normal density p{r) is shown in Fig. 5, the probability density P{r) — 47rr^p(r) 
on a log scale in is shown Fig. 6, and the probability density on a linear scale is shown Fig. 
7. In all figures the neutron density is shown by the dashed line, the core proton density 
with crosses, and the valence proton density with a sohd fine. Note in Figs. 6—7 that the 
areas are equal to three, four and one, respectively. The valence proton clearly has a large 
extension, but whether or not it constitutes a "halo" or a "skin" is a question of semantics. 

The valence proton (vp), core proton (cp) and neutron (n) rms radii are: 




(21) 




(22) 



and 




(23) 



The total proton rms radius is given by 




(24) 



The rms charge radius includes the rms radius of 0.80 fm for the proton: 




(25) 



and the matter radius is given by: 




(26) 
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Our results can be compared to those of other theoretical calculations. The results 



of \J{r'^)vp = 3.75 fm obtained by Riisager and Jensen is much smaller than ours but they 
use an arbitrary potential shape. The reason for their small 5*17 is obvious from Fig. 3. The 
calculations presented in Ref 10 for the Q moment appear to be very close to our results. 

The interaction cross sections for and ^Li on ^^C have been calculated with the 
method of Ref 34. The results with the finite-range interaction are cr=843 mb for and 
(7=820 mb for ^Li. These can be compared to the experimental values of (7=784(14) mb^^ 
for ^B and (t= 768(9) mb^^ for ^Li. The agreement for the magnitudes is as good as can 
be expected from the uncertainties in the calculation. However, the ^B/^Li ratio (in which 
some of the reaction uncertainties may cancell) is in excellent agreement between theory and 
experiment. The effects of the proton "halo" in ^B and the neutron "halo" in ^Li are not 
large compared with the classic cases ^'''^^ of ^^Li, ^^Be and ^''Be. 

Recent radiaoactive beam experiments ^^'^^ have looked at the momentum distribution 
of the ^Be fragments which result from the break-up of a beam ^B on various targets. Prom 
these experiments one expects to measure the momentum distribution of the most loosely 
bound protons. The longitudinal momemtum P{kz) obtained for the valence proton is shown 
in Fig. 8. This is obtained by the Fourier transform of the spacial wave function: 





(27) 



where ^'(r) = ip{r)Y(_{f) '^{k) = ip{k)Yi{k), and where the radial momentum distribution is 



given by: 




(28) 



The logitudinal momentum distribution is given by 




(29) 



where k"^ = k^ + kl- 
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The calculated momentum distribution has a width of about 150 MeV/c compared 
to the experimental value of 81±6 MeV/c. Given the good agreement generally found for 
calculated and obseved neutron halos^^ this disagreement is puzzling. There is some dis- 
cussion in the literature about ways to improve the above calculation to take into account 
the peripheral nature of the reaction. ^^'^'^ In the peripheral direct reaction model one puts 
in an additional cut-off in Eq. (29) to exclude the interior part of the radial distribution 
which presumably does not contribute because the cross section coming from that part is 
dominated by a more violent reaction where the core (^Be in this case) is destroyed. We 
have phenomenologically modeled this effect by putting a Fermi shaped cut-off factor in Eq. 
(28) which has the effect of excluding the interior out to a radius Rcut and with a diffuseness 
(icut- We take acut=0.65 fm and vary Rcut to get about the observed momentum distribtion. 
This requires Rcut=5 fin and the results are shown by the dashed line in Fig. 8. The cut-off 
results in a reduction of P{kz) at small momenta by a factor of 2.5 and we have renormalized 
the cut-off distribution by this factor in order to show the change in width. It is already 
known that the cut-off factor does not have much effect on the neutron halo momentum 
distributions,^^ and we have demonstated that even a value as large as Rcut = 5 fm has 
little effect on the width of the ^^Be neutron halo momentum distribution. We do not know 
why we should need Rcut=^ but we note that this corresponds to the point in Fig. 5—7 
where the valence proton density falls below the core density. More work needs to be done 
to understand these results. 

Schwab et al. "^^ present an RPA calculation which goes beyond the Op shell and which 
agrees with the shape of the observed momentum distribution. However, the shape of the 
wave function beyond about 6 fm as shown in Fig. 4 in their paper appears unrealistic to 
us. The shape beyond about 6 fm is entirely determined by the Coulomb and centrifugal 
barriers, and their shape is differs from this expectation. Our own calculation can of course 
be criticized for staying within the Op shell. However, a very recent "no-core" calculation^^ 
along the lines of those given in Ref 42 which takes into account the lowest six major shells 
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(21 shell-model orbitals) and up to Ahcu in excitation gives spectroscopic factors which are 
close to the present Op shell results. 

VII. CONCLUSIONS 

Our calculated value of 5'i7(20keV) = 26.5 ± 2.0 eV-barns is agreement with the 

higher values of 25-27 eV-barns inferred from the (p,7) data of Parker^ and Kavanagh et 

al., ^ and is higher than the value obtained from the weighted average of all experimental data, 

Sn — 22.2 ± 2.3 eV-barns, which is the value currently adopted in most solar models. 

As far as other theoretical predictions are concerned, our current result for 5*17 is roughly 

near the highest of these (~30 eV-barns) '^^ and and much larger than the smallest (~17 
eV-barns). 24.46 

We note, that in contrast to the common behef (e.g. Ref 46), a small value of 5'i7 
does not make the solar neutrino problem less severe. If one takes standard nuclear and solar 
physics, and standard neutrino properties, then the best fit of the neutrino fluxes indicates 
a suppression in both the ^Be (^7) and the (^g) neutrino fluxes, but the suppression is 
much stronger in ^7. However, a smaller 5'i7 value alone would make the predicted 07/08 
ratio larger, and hence would exagerage the solar neutrino problem. 
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Figure Captions 

Fig. 1: The valence proton rms radius as a function of Aq. The three hues join the points 
obtained for different value of the diffuseness a for R=2.0 fm (filled circles), i?=2.4 fm (open 
circles) and R—2.S fm (squares). The results for the specific potentials in Table II are shown 
by the labeled crosses. 

Fig. 2: The valence proton density at r=10 fm as a function of Aq (see caption to Fig. 1). 

Fig. 3: The valence proton density at r=10 fm as a function of the valence proton rms radius 
(see caption to Fig. 1). 

Fig. 4: The scattering wave function for ''Li(n,7) obtained with the standard potential (BI) 
of Barker (sohd hne) together with the results obtained with the cluster model (dashed 
line) . 

Fig. 5: The radial density profile for ^B for the neutrons (dashed hne), the core protons 
(crosses) and the valence proton (solid line). 

Fig. 6: The radial probability distribution for ^B on a log scale (see caption to Fig. 5). 

Fig. 7: The radial probability distribution for "^B on a linear scale (see caption to Fig. 5). 

Fig. 8: The momentum distribution for the ®B valence proton. The solid line corresponds 
to the full radial wave function and the dashed line corresponds to the radial wave function 
cut-off at r=5 fm. 
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TABLES 



TABLE I. 
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TABLE II. Woods-Saxon potential parameters from the present analysis (sets A, B and C) 
compared to those used by Tombrello^^ and Barker I. Barker II and III correspond to those 
values Barker needed to reproduce the ^Li(n,7) cross section. 

Set R (fm) a (fm) 

A 2.8 0.51 

B 2.4 0.69 

C 2.0 0.81 

Tombrello (T) 2.95 0.52 

Barker BI 2.39 0.65 

Barker BII 0.53 0.65 

Barker Bill 2.39 0.27 
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TABLE III. Spectroscopic amplitudes 9j{Ji, Jf = 3/2) for to ^Be from the CKl/^ Kumar 

and PTBME^^ interactions. The order of the states for a given JJ is indicated by n/. 

JJ,nf j CKI Kumar PTBME 

2+,l 3/2 0.988 0.966 0.986 

1/2 -0.237 -0.259 -0.253 

1+,1 3/2 0.567 0.606 0.552 

1/2 -0.352 -0.244 -0.342 

3+,l 3/2 0.581 0.555 0.565 

l+,2 3/2 0.617 0.574 0.525 

1/2 0.840 0.861 0.859 
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TABLE IV. Spectroscopic factors S = S1/2 + S'3/2 for ^Li and ^B. 
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